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ABSTRACT 

Pulsars have been recognized as normal neutron stars, but sometimes argued as quark stars. 
Sub-millisecond pulsars, if detected, would play an essential and important role in distin- 
guishing quark stars from neutron stars. We focus on the formation of such sub-millisecond 
pulsars in this paper. A new approach to form a sub-millisecond pulsar (quark star) via ac- 
cretion induced collapse (AIC) of a white dwarf is investigated here. Under this AIC process, 
we found that: (1) almost all the newborn quark stars could have an initial spin period of 
~ 0.1 ms; (2) the nascent quark stars (even with a low mass) have sufficiently high spin- 
down luminosity and satisfy the conditions for pair production and sparking process to be as 
sub-millisecond radio pulsars; (3) in most cases, the timescales of newborn quark stars in the 
phase of spin period < 1 (or < 0.5) ms can be long enough to be detected. 

As a comparison, an accretion spin-up process (for both neutron and quark stars) is also 
investigated. It is found that, quark stars formed through AIC process can have shorter periods 
(^ 0.5 ms); while the periods of neutron stars formed in accretion spin-up process must be 
longer than 0.5ms. Thus if a pulsar with a period less than 0.5 ms can be identified in the 
future, it should be a quark star. 
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1 INTRODUCTION 

Though it has been more than 40 years since the discovery of ra- 
dio pulsars, their real nature is still not yet clear because of the 
uncertainty about cold matter at supranuclear density. Both neu- 
tron matter and quark matter are two conjectured states for such 
compact objects. The objects with the former are called neutron 
stars, and with the latter are quark stars. It is an astrophysical chal- 
lenge to observationally distinguish real quark stars from neutron 
stars (see reviews by, e.g., Madsen 1999; Glendenning 2000; Lat- 
timer & Prakash 2001; Kapoor & Shukre 2001; Weber 2005; Xu 
2008). The most obvious discrepancy could be the minimal spin 
period of these two distinct objects. The minimal periods of these 
two kinds of objects are related to their formation process. How 
fast a neutron star or a quark star can rotate during the recycling 
process in low mass X-ray binary (LMXB) has been considered 
by several authors (Bulik, Gondek-Rosinska & Kluzniak 1999; 
Blaschke et al, 2002; Zdunik, Haensel & Gourgoulhon 2002; Xu 
2005; Arras 2005). Friedman, Parker & Ipser (1984) have found 
that neutron stars with the softest equation of state can rotate as 
fast as 0.4 ms. The smallest spin period for neutron stars computed 
by Cook, Shapiro & Teukolsky (1994) is about 0.6 ms. Frieman 
& Olinto (1989) have showed that the maximum rotation rate of 
secularly stable quark stars may be less than 0.5 ms. Burderi & 
D'Amico (1997) have discussed a possible evolutionary scenario 
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resulting in a sub-millisecond pulsar and the possibility of detect- 
ing a sub-millisecond pulsar with a fine-tuned pulsar-search survey. 
Gourgoulhon et al. (1999) have investigated the maximally rotating 
configurations of quark stars and showed that the minimal spin pe- 
riod was between 0.513 ms and 0.640 ms. Burderi et al. (1999) 
have predicted that there might exist a yet undetected population 
of massive sub-millisecond neutron stars, and the discovery of a 
sub-millisecond neutron star would imply a lower limit for its mass 
of about 1.7 Mq. A detailed investigation about spin-up of neu- 
tron stars to sub-millisecond period, including a complete statistical 
analysis of the ratio with respect to normal millisecond pulsars, was 
performed by Possenti et al. (1999). The minimal recycled period 
was found to be 0.7 ms. Gondek-Rosinska et al. (2001) have found 
that the shortest spin period is approximately 0.6 ms through the 
maximum orbital frequency of accreting quark stars. Huang & Wu 
(2003) have found that the initial periods of pulsars are in the range 
of 0.6 ~ 2.6 ms using the proper motion data. Zheng et al. (2006) 
have showed that hybrid stars instead of neutron or quark stars may 
lead to sub-millisecond pulsars. Haensel, Zdunik & Bejger (2008) 
have discussed the compact stars' equation of state (EOS) and the 
spin-up to sub-millisecond period, via mass accretion from a disk 
in a low-mass X-ray binary. 

There have been many observational attempts in searching 
sub-millisecond pulsars. A possible discovery of a 0.5 millisecond 
pulsar in Supernova 1987A is not held true in the follow-up ob- 
servations (Sasseen 1990; Percival et al. 1995). Bell et al. (1995) 
reported on optical observation of the low mass binary millisecond 
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pulsar system PSR J0034-0534, and they used white dwarf cooling 
models to speculate that, the limit magnitude of the J0034-0534's 
companion suggested that this millisecond pulsar's initial spin pe- 
riod was as short as 0.6 ms. As addressed by D'Amico & Burderi 
( 1999), in particular the detection of a pulsar with a spin period well 
below 1 ms could put severe constraints on the neutron star struc- 
ture and the absolute ground state for the baryon matter in nature. 
They designed an experiment to find sub-millisecond pulsars with 
Italian Northern Cross radio telescope near Bologna. Edwards, van 
Straten & Bailes (2001) have found none of sub-millisecond pul- 
sars in a search of 19 globular clusters using the Parks 64 m Radio 
telescope at 660 MHz with a time resolution of 25.6 (is. Han et al. 
(2004) did not find any sub-millisecond pulsars from highly polar- 
ized radio source of NVSS (NRAO VLA Sky Survey). Kaaret et al. 
(2007) have found oscillations at a frequency of 1 122 Hz in an X- 
ray burst from a transient source XTE J 1739-285 which may con- 
tain the fastest rotating neutron star so far. Significant difficulties 
do exist in current radio surveys for binary sub-millisecond pulsars 
due to strong Doppler modulation and computational limitations 
(Burderi et al. 2001). 

How do sub-millisecond pulsars form? This is still an open 
question which we will explore in this paper. Previously, discus- 
sions were concentrated on the formation of neutron stars or quark 
stars spun up via accretion in binaries. We have considered a new 
approach to create a sub-millisecond pulsar (quark star) with super- 
Keplerian spin via accretion induced collapse (AIC) of a mas- 
sive white dwarf (WD). The initial spin of the newborn quark star 
could be super-Keplerian, and it can have a long lifetime in sub- 
millisecond phase and produce enough strong radio emission to be 
detected. 

In §2, we discuss low mass quark stars formed from AIC of 
WDs, which can have an minimal initial period of sub-millisecond. 
In §3, the radiation parameters and the conditions for pair produc- 
tion are estimated in order to investigate whether the AlC-induced 
quark stars could be pulsars or not. The lifetimes of sub-millisecond 
pulsars are also estimated and the possibility of detection is dis- 
cussed. The spin-down evolution diagrams of a newborn quark star 
and neutron star are also plotted. In §4, as a comparison, the sub- 
millisecond pulsars formed through accretion acceleration (spin- 
up) in binary systems are also considered. In §5, conclusions and 
discussions are presented. 



2 SUB-MILLISECOND QUARK STARS FORMED 
THROUGH AIC OF WHITE DWARFS 

Neutron star's formation from AIC of a massive white dwarf is 
widely discussed by many authors (Nomoto et al. 1979; Nomoto 
& Kondo 1991; van Paradijs et al. 1997; Fryer et al. 1999; Bravo 
& Garcfa-Senz 1999; Dessart et al. 2006). Recently, it is pointed 
out that Galactic core-collapse supernova rate cannot sustain all the 
separate neutron star populations (Keane & Kramer 2008), which 
implies other mechanisms for forming neutron stars. AIC of a mas- 
sive WD can be an important mechanism for pulsar formation, even 
for isolated pulsars if the binary systems are destroyed due to strong 
kicks. We now discuss the possibility for a low mass quark star 
formed from AIC of a WD. In a binary system, when the WD 
has accreted enough matter from the companion so that its mass 
reaches the Chandrasekhar limit, the process of electron capture 
may induce gravitational collapse. The detonation waves burn nu- 
clear matter into strange quark matter which spread out from the 
inner core of the WD (Lugones, Benvenuto & Vucetich 1994). A 
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Figure 1. The relation of mass and radius of WDs. The red line is theoreti- 
cal line. The blue triangles and circles are observed WDs' data which were 
taken from Table 1 of Nalezyty & Madej (2004) and Table 3 & 5 of Proven- 
cal et al. (1998), respectively. Among these data, the WD RE J0317-853 is 
the most massive WD, whose mass and radius are 1.34Mq and 2400 km 
respectively. The square is the point of a WD with the Chandrasekhar mass 
limit. 

boundary of strange quark matter and nuclear matter will be found 
at the radius where the detonation waves stop when nuclear mat- 
ter density drops below a critical value. A similar process was also 
discussed and calculated by Chen, Yu & Xu (2007). The size of 
the inner collapsed core may depend on the chemical composition 
and accretion history of the WD (Nomoto & Kondo 1991). Conse- 
quently, quark stars with different masses could be formed. 

Both rigidly and differentially rotating WDs are taken into ac- 
count. As a first step, we assume that both the collapsed WD and 
the newborn quark star have rigidly rotating configurations for sim- 
plicity. The WDs, progenitors of these quark stars, could have a uni- 
formly rotating configuration due to the effects of crystallization, 
as well as an increase of central density may lead to catastrophic 
evolution (supernova) (Koester 1974). With these assumptions, a 
model of sub-millisecond pulsars' formation is given below. The 
initial spin period of AlC-produced quark stars can be estimated as 
follows. We assume that the mass (M*) of the nascent quark star 
ranges from 1O _3 M0 to IM@, and the white dwarf rotates rigidly 
at the Kepler period (Pk) just before collapsing. The quark star's 
rest mass (M*) is approximately equal to the mass (m CO re) of the 
inner collapsed core of the white dwarf. If the angular momentum 
is conserved during AIC, the newborn quark star can rotate at a 
much shorter period, P q , then 

/ *L-i*L (i) 

J core p — Jq . \l) 
±K -i i q 

This is to say, 

A = -^Pk, (2) 

J core 

where 7 q is the quark star's moment of inertia, and J C oro is the 
moment of inertia of WD's inner collapsed core, which can be well 
approximated by 
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(3) 



The mass and radius of a low mass (A/* ^ IMq) quark star could 
be approximately related by M* = (4/3)7r(4,3)P 3 (Alcock, Farhi 
& Olinto 1986) in the bag model. We have an approximate formula 
for the fast rotating quark star's moment of inertia 
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where z-axis is spin axis; x is integral variable of each disc perpen- 
dicular to the spin axis; c is the speed of light; R is the quark star's 
radius; the bag constant f3 of quark stars is (60-110) MeV fm -3 , 
i.e. (1.07-1.96) x 10 14 g cnT 3 , /3 14 in units of 10 14 g cm~ 3 . 

For WD, we made a code, using both non-relativistic hydro- 
static equilibrium equation 



dp 
dr 



Gm(r)p(r) 



(5) 



and general equation of state (EOS) for a completely degenerate 
fermi gas 



37T 2 ft 3 



2 4 

c p 



y/c 2 p 2 + m 2 cf 
= 1.42 x 10 25 (j>(x) dyn cm~ 2 , 



:dp 



(6) 



where x = pf/itic, A e = fi/(mc) the electron's Compton wave- 
length, Pf the fermi momentum, 

4>{x) = {8tv 2 )- 1 {x(1 + x 2 ) 1/2 (2x 2 /3-1) 
+ \n[x + (1 + x 2 ) 1/2 ]} 

to calculate the mass (m CO ro) and moment of inertia (/core) of the 
collapsed core of a massive WD, where p, p, G, h are pressure, 
mass density, the gravitational constant and the Planck constant, 
respectively. 

Using Eqs. (5) and (6), one can make numerical calculation to 
get the WD's theoretical relation of mass and radius (the red line 
in Figure 1). For comparison, one can see a figure on lineQ. Before 
collapsing, the mass of a WD is close to the Chandrasekhar mass 
limit, as high as Mwd = 1.4Mq (Shu 1982), the corresponding 
radius is much smaller, such as Pwd = 410 km. We could nu- 
merically obtain the initial period, P q , of nascent quark stars with 
different mass via Eq. (2), and find that almost all the values of P q 
are around ~ 0.1 ms (See Table 1) if the WD rotates rigidly at an 
almost Kepler period due to accretion (or spin-up) in a binary just 
before collapsing. The newborn quark stars' surface spin velocities 
are well above the Kepler velocities, we regard this as "the super- 
Keplerian case". 

WDs may be rotating differentially. The detailed calculations 
are given in Appendix A. Therefore, as a follow-up second step, we 
also use Eqs. (2), (5), (6) and (A2) to calculate the initial spin period 
of the nascent quark stars in the differentially rotating WD model, 
taking the free parameter a = 0.5. The results are also shown in 
Table 1. 

1 |http : //cococubed. asu ■ edu/co de-pages/coldwd . shtml| 



A newborn quark star could certainly rotate differentially, and 
may be relaxed to become a rigidly rotating configuration finally. 
However, the timescale of the relaxa tion depen ds on the viscosity 
and the state of cold quark matter dXiil2009l) . Nevertheless, the 
newborn quark star's relaxation (from differentially rotating con- 
figuration to rigidly rotating configuration) may be due to fast so- 
lidification after birth. A calcu lation shows that t he solidification 
timescale is only 10 3 - 10 6 s jXu & Liang||2009h . Therefore, the 
relaxation timescale could be much shorter than the lifetime of pul- 
sars within sub-millisecond periods (See the following section 3.3). 

The WD RE J0317-853 has the highest observed mass (1.34 
M© close to the Chandrasekhar limit) with radius of 2400 km 
(Nalezyty & Madej 2004). If a WD like RE J03 17-853 could be in 
a binary and accreted enough materials to the Chandrasekhar limit, 
then it may collapse. Therefore, under this assumption, we also cal- 
culated the initial spin periods P q and Pdif of a nascent quark star. 
The calculated results are listed in Table 1 . It is found that, even if 
a WD has a larger radius such as 2400 km, it can also collapse to a 
sub-millisecond quark star for either rigidly or differentially rotat- 
ing WD models. In the differentially rotating WD model, it tends 
to give a rigidly rotating configurations in the limit of large values 
of a, Pdif increases as the parameter a increases. The conclusions 
from the rigid rotation model are valid even if differential rotation 
is included. 



Can a quark star survive even if it rotates at such a high fre- 
quency (~ 10 4 Hz)? Will it be torn apart by the centrifugal force? 
There are quite distinguishing characteristics between neutron stars 
and quark stars. A low mass quark star is possible to spin at a super- 
Keplerian frequency because it is self-bound by strong interaction. 
On one hand, as noted by Qiu & Xu (2006), astrophysical quark 
matter splitting could be color-charged if color confinement cannot 
be held exactly because of causality. On the other hand, however, 
rapidly spinning quark matter could hardly split if color confine- 
ment is held exactly. In addition, the recently discovered nature of 
strongly coupled quark gluon plasma (sQGP) as realized at Rela- 
tivistic Heavy Ion Collider (RHIC) experiment (e.g., Shuryak 2006) 
may also prevent a super-Keplerian quark star to split. 

The short spin period above is not surprising, and could be ver- 
ified for a simplified special case, if both quark star's density (=4/3) 
and white dwarf's density (= p c ) are uniform. Using Eq. (2) and 
the mass-radius relation, we can find the initial period of the quark 
star to be P q = O c /4/3) 2/3 P WD ~4x 10- 3 (pu//3i 4 ) 2/3 PwD 
(with Pwd the spin period of white dwarf, pw = p c /10 11 g cm 3 , 
/3i4 = /3/10 14 g cm 3 ), which depends only on the densities of WD 
and quark star. 



If the WD has not been spun up fully to the Kepler period, i.e., 
the WD rotates at a sub-Keplerian period (e.g., several times of Pk) 
before AIC, can the initial period of a newborn quark star formed 
from such a WD be of sub-millisecond? We investigated the case 
of a massive WD (1.4Mq, 410 km) rotating at a period Pwd = 
5Pk ~ 600 ms. The initial spin period of quark stars with different 
mass are as follows: P q ~ 0.11 ms for a quark star with mass 
of 0.001M Q ; ~ 0.24ms for O.OLM ; 0.35ms for O.1M and ~ 
0.36 ms for \Mq. The spin-down feature of such a newborn quark 
star depends on its gravitational wave radiation and magnetodipole 
radiation (see details in §3). 
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Table 1. The minimal initial period (P q ) and lifetimes (r) due to GW and EM radiation in the phase of sub-millisecond 
period for quark stars with different masses in the super-Keplerian case. P q and P^if are calculated via angular momentum 
conservation using rigidly and differentially rotating WD model with central density of 10 11 gcm — 3 . P q and P di[ are 
similarly calculated but using a WD like RE J0317-853 with mass of IAMq and radius of 2400 km. t\ is quark stars' 
lifetime in the phase of < 1 ms, while T2 is the timescale in the phase of < 0.5 ms. P q is also used in Table 2 & 3 and Figure 
2. (3 is the bag constant, e e is the gravitational ellipticity. 



Mass 

{M e ) 


Radius (km) 
/3 = 60 MeV fm~ 3 


Pq(ms) 
(3 = 60 MeV fm~ 3 


P dif (ms) 
a = 0.5 


Pq(ms) 
(3 = 60 MeV fm" 3 


Piif(ms) 
a = 0.5 


0.001 


1.04 


0.0699 


0.0261 


0.0481 


0.0252 


0.01 


2.24 


0.0751 


0.0472 


0.0512 


0.0470 


0.1 


4.81 


0.104 


0.101 


0.102 


0.101 


1 


10.37 


0.221 


0.218 


0.218 


0.218 



Mass(M Q ) 


P q (ms) 


n(yr) 

£e = 10" 6 Be = 10" 9 


T2(yr) 

e e = 10" 6 e e = 10" 9 


0.001 


0.0699 


3.4 x 10 7 


4.5 x 10 10 


2.1 x 10 e 


1.1 x 10 10 


0.01 


0.0751 


7.3 x 10 s 


2.0 x 10 10 


4.5 x 10 4 


4.3 x 10 9 


0.1 


0.104 


1.6 x 10 4 


5.4 x 10 9 


9.0 x 10 2 


6.5 x 10 8 


1 


0.221 


3.4 x 10 2 


3.1 x 10 s 


2.2X10 1 


2.0 x 10 7 



3 RADIATION OF SUB-MILLISECOND QUARK STARS 
WITH LOW MASSES 

The mass of most sub-millisecond quark stars formed from WD's 
AIC is so low, can the quark star produce radiation luminous 
enough to be observed like millisecond pulsars? This is related 
to two aspects. First of all, is the rotational energy loss rate high 
enough to power the electromagnetic radiation as normal neutron 
stars? Secondly, is the potential drop in the inner gap high enough 
for pair production and sparking to take place in the inner gap? 
These are necessary conditions for radio emission of pulsars. 

3.1 The spin-down power of sub-millisecond pulsars 

Normal radio pulsars are rotation-powered, and the radiation en- 
ergy is coming from the rotational energy loss. Here we neglect 
gravitational wave radiation first, then the rate E Iot , is 

■ 8TY i R 6 B 2 p- 4 

^ = 3? ■ (7) 

Comparing the rotational energy loss rate (E ro t,q) of quark stars 
with normal neutron stars' (-Erot.Ns), one can have 



-E^rot.q/ E?, 



RqBqPq 

r>6 d2 p-4 



(8) 



If we take normal parameters, such as the surface magnetic field 
of polar cap B q = 10 s G, Bns = 10 12 G, the rotational period 



0.1 ms and Pns = 1 s, the result is E I0 t,q/E- 



rot,NS 



10" 



even for a quark star with a mass of 0.001 Mq . This means that the 
quark stars have enough rotational energy to radiate, hundred times 
than normal pulsars, even if the mass is so low. 

3.2 Particle acceleration for sub -millisecond pulsars 

In most radio emission models of pulsars, such as RS model (Rud- 
erman & Sutherland 1975, hereafter RS75), inverse Compton scat- 
tering (ICS) model (Qiao & Lin 1998), the multi-ring sparking 
model (Gil & Sendyk 2000), the annular gap model (Qiao et al. 
2004) and so on, the potential drop in the inner gap must be high 
enough so that the pair production condition can be satisfied. 



In the inner vacuum gap model, there is strong electric field 
parallel to the magnetic field lines due to the homopolar generator 
effect. The particles produced through 7 — B process in the gap can 
be accelerated to ultra-relativistic energy (i.e., the lorentz factor can 
be 10 6 for normal pulsars). The potential across the gap is (RS75) 



A V 



h , 

c 



(9) 



where SI is the angular frequency of the pulsar; h is the gap height; 
B and c represent the magnetic field at the surface of the neutron 
star and the speed of light, respectively. As h increases and ap- 
proaches r p , the potential drop along a field line traversing the gap 
can not be expressed by Eq. (9) above. In this case the potential can 
reach a maximum value 



A V m 



OB 2 

2c Vl " 



(10) 



where r p is the radius of the polar cap. 

Let us make an estimate about the quark star's potential drop 
AV q in the polar gap region 

QB a 



AV q = 



2c rp ' q ' 



(11) 



where D. = 2n/P q , r p , q = i? q (27r_R q /cP q ) 1//2 . For normal neu- 
tron stars, AV can be obtained by just changing the subscript q to 
NS. Thus 



AV q B^P- 



AVns BnsRI^Pk 2 ' 



(12) 

l NS J NS 

As one can take R q = 1 km for a quark star with the mass of 
0.001M©, B q = 10 s G, B NS = 10 12 G, R NS = 10 km, P q = 
0.1 ms and P NS = 1 s, we find that AVq/ A V"ns = 10. This 
means that the quark stars can have enough potential drops in the 
polar cap regions. 

In the inner gap model, 7 — B process plays a very impor- 
tant role, two conditions should be satisfied at the same time for 
pair production: (1) to produce high energy 7-ray photons, a strong 
enough potential drop should be reached; (2) for pair production, 
the energy component of 7-ray photons perpendicular to the mag- 
netic field must satisfy _E 7 .x ^ 2m e c 2 (Zhang & Qiao 1998). 

Particles produced in the gap can be accelerated by the electric 
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Table 2. Gap parameters estimated for sub-millisecond quark stars. E TO t is the spin-down luminosity; /icr ' s me curvature 
radiation (CR) gap height; AVqr is the potential drop of CR gap; h IaB is the height of resonant ICS gap; AV r es is the 
potential drop of resonant ICS gap; h tri is the thermal ICS gap height; A V t h is the potential drop of thermal ICS gap. 



M(M Q ) 


I?rot (erg s 1 ) 


/iCR(cm) 


AVcr(V) 


/l res (cm) 


AV res (V) 


h th (cm) 


AV th (V) 


0.001 


4.99 x 10 36 


1.16 x 10 4 


2.84 x 10 10 


3.13 x 10 5 


2.05 x 10 13 


1.18 x 10 3 


2.93 x 10 8 


0.01 


3.75 x 10 38 


1.34 x 10 4 


3.76 x 10 10 


3.64 x 10 5 


2.78 x 10 13 


1.31 x 10 3 


3.62 x 10 s 


0.1 


1.02 x 10 40 


1.72 x 10 4 


6.19 x 10 10 


4.61 x 10 5 


4.46 x 10 13 


1.62 x 10 3 


5.47 x 10 s 


1 


5.00 x 10 40 


2.65 x 10 4 


1.47 x 10 11 


6.74 x 10 5 


9.52 x 10 13 


2.36 x 10 3 


1.16 x 10 9 



field in the gap and the Lorentz factor of the particles can be written 

as 

eAV 



7 



.12 ' 



(13) 



where 7 is the Lorentz factor of the particles accelerated by the 
potential A V, m e the mass of an electron or positron, e the charge 
of an electron. 

In 7 — B process, the conditions for pair production are that 
the mean free path of 7-ray photon in strong magnetic field is equal 
to the gap heights, / ~ h. The mean free path of 7-ray photon is 
given by (Erber 1966) 



(14) 



, _ 4.4 g B c j4_ 

eyhcm e cB ± em 3 X 

where B c = 4.414 x 10 13 G is the critical magnetic field, h the 
Planck's constant, 



X 



B 



E-y B\ 



- sin v — 

2m e c 2 B c 



(15) 



2m e c 2 B c ' 

and B± is the magnetic field perpendicular to the moving direction 
of 7 photons, which can be expressed as (RS75) 

B ± ^-B^-B. (16) 
P P 

Here I w h is the condition for sparks (pair production) to take 
place, p is curvature radius of the magnetic field lines. For dipole 
magnetic configuration, it is (Zhang et al. 1997a) 



(17) 



where A is a parameter to show the field lines, A = 1 corresponding 
to the last opening field line. Gamma-ray energy from the curvature 
radiation process can be written as 



£ /ry £j-> 



37 3 c 
2p ■ 



(18) 



We estimated the gap heights based on Zhang, Qiao & Han 
(1997b), i.e. 

,6 ; 



% , „o r>3/7 D-4/7 2/7 

/icr ~ 10 P ' B g p 6 ' cm 



(19) 



When the relevant parameters used are B = 10 G , P = P q 
and assuming a dipole magnetic configuration, for any mass quark 
stars, one can estimate the gap height from curvature radiation (CR) 
h cr « 10 4 cm = 100 m. This means that even if without multipo- 
lar magnetic field assumption, the quark star can still work well for 
the CR pair production. 

There are three gap modes for pair production, i.e. resonant 
ICS mode, thermal-peak ICS mode and CR mode (Zhang et al. 
1997a). Each mode has relevant gap parameters including gap po- 
tential drop AV and the mean free path I of 7 — B process. For nor- 
mal neutron stars, one needs the assumption of a multipolar mag- 
netic field, p = 10 6 cm, as RS75; but for 0.1 ms low mass quark 



stars, the dipole curvature radius is about 10 cm. We estimated gap 
heights and other parameters based on the work of Zhang, Qiao & 
Han (1997b), as shown in Table 2. 

One can see from Table 2 that when the high energy gamma- 
ray photons come from resonant photon production, the height of 
the gap is larger. For the thermal-peak ICS mode, it is one order of 
magnitude lower than the CR mode, and two order of magnitude 
lower than resonant ICS mode. This means that in most cases, the 
thermal-peak ICS induced pair production is dominated in the gap. 

The newborn sub-millisecond quark stars have enough spin- 
down luminosities and gap potential drops (see Table 2), so that 
they may emit radio or 7-ray photons with sufficient luminosities, 
which can be detected by new facilities, e.g., FAST and Fermi (for- 
merly GLAST). 

3.3 Lifetimes of the sub-millisecond pulsars in the phase of a 
short spin period 

Sub-millisecond pulsars may be very rare, or the timescale for such 
a pulsar to stay in the short period phase (< 1 ms) may not be 
long enough due to magnetodipole (EM) radiation and gravitational 
wave (hereafter GW) radiation (Andersson 2003). The lowest or- 
der GW radiation is bar-mode, which is due to non-axisymmetric 
quadrupole moment. Here we consider GW radiation on the bar 
mode which exerts a larger braking torque with braking index 
n ~ 5 than magnetodipole radiation (n = 3). The rotation fre- 
quency drops quickly due to GW radiation and EM radiation: 



5 C s ' 6c 3 ' (20) 

where c is the speed of light, e e = Aa/a is the gravitational ellip- 
ticity (equatorial ellipticity), Aa is the difference in equatorial radii 
and a is the mean equatorial radii. 

To simplify Eq. (20), we introduce the notation A = 
32G/£e/(5c 5 ) and D = B^R 6 /(61c 3 ), and integrate the equation 
in the angular velocity's domain [fij = 27r/Pj,f2o = 27r/0.001], 
then 



= —( — 



fi?' 



A 
2D 2 



In 



A 
D 



1 , A 



(21) 



An accurate ellipticity of quark stars is unfortunately uncer- 
tain. Nevertheless, let's estimate the e e to calculate the timescales 
in the sub-millisecond period phase for GW and EM radiations. 
Cutler & Thorne (2002) suggested e e = (/ - I )/I ^ 10" 6 . 
Regimbau & de Freitas Pacheco (2003) found from their simu- 
lations that e e — 10~ 6 is the critical value to have an at least 
one detection with interferometers of the first generation (LIGO or 
VIRGO). It was shown that direct upper limit was e e ~ 1.8 X 10~ 4 
on GW emission from the Crab pulsar using data from the first 
9 months of the fifth science run of LIGO (Abbott et al. 2008). 
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In addition, Owen (2005) showed that the maximum ellipticity of 
solid quark stars was £ e ,max = 6 x 10~ 4 . From the on-line cata- 
logue hosted by the ATNF0, the seventh fastest rotating millisec- 
ond pulsar is PSR J0034-0534, which has very low period deriva- 
tive P ~ 4.96 x 10~ 1 ss _1 . We thus use such a low P and 
Eq. (20) to constrain the lower limit of the sub-millisecond pul- 
sars' ellipticity, which is e e ,mm ~ 10~ 9 if the stellar mass is one 
order of one Solar mass. For quark stars, in order to facilitate to 
compare with the neutron stars' lifetime (r) in the phase of sub- 
millisecond period, we use mean equatorial ellipticities e e — 10~ 6 
and £ e = 10 to calculate r for both quark stars and neutron stars 
through Eq. (21). 

In the case of e e = 10" 6 , if we make the hypothesis that the 
rotational energy is lost because of EM radiation, then one can eas- 
ily derive t em = l/(2£>)(l/fi§ - ~ 5.9 x 10 9 yr for a 
typical compact star with Bo ~ 10 8 G and M = Mq . While, if we 
suppose that the rotational energy is lost due to GW radiation, then 
row = l/(4A)(l/fi - 1/fi?) ~ 10 2 yr for a typical compact 
star. The energy loss rate of GW & EM radiation in the phase of 
sub-millisecond period, for a typical compact star which has a low 
magnetic field (10 8 -10 9 G) either from AIC (Xu 2005) or spun up, 
are E G w = 32GI 2 e 2 f2 6 /(5c 5 ) = 7.0 x 10 41 P^ a 6 ergs" 1 and 
Eem = B 2 i? 6 fi 4 /(6c 3 ) = 9.6 x 10 34 P~ 4 B 2 i?| ergs" 1 , respec- 
tively. Even if a quark star with 1 Mq formed from WD's AIC has 
a high magnetic field such as 10 12 G, the lifetime r in the phase of 
sub-millisecond is 37 years, in comparison with r = 336 years for 
Bo = 10 s G. Then the EM energy loss is similar to the GW energy 
loss and becomes very important for Bo = 10 12 G. For Bo ranges 
from 10 8 G to 10 11 G, one always has Bgw 3> Eem for compact 
stars with short spin periods (< 1 ms). Therefore, in the case of 
larger ellipticity (e. g., e e = 10 -6 ), it is clear that GW radiation 
dominates the energy loss in the phase of short period for either 
recycled or AIC's compact stars with low magnetic field. The cor- 
responding lifetime is shorter for a compact star with higher mass 
(~ Mq), but longer for a star with lower mass (~ O.OOIMq). 
However, if the ellipticity is lower, such as e e = 10~ 9 , EM radi- 
ation dominates the rotational energy loss. The corresponding life- 
time of a quark star (even with a high mass ~ Mq) is long enough 
for us to detect. Figure 2 shows the relation of lifetime (in the phase 
of < 1 ms) and gravitation ellipticity e c for quark stars. 

In the super-Keplerian case, the timescales in the phase of 
<0.5 ms for quark stars with different mass are also calculated, 
and listed in Table 1 (See r-£). For a high-mass quark star with 
larger ellipticity, the timescale is too small for real detection; but 
the timescale is > 10 4 yr for a low mass quark star. Therefore, low 
mass quark stars with e e ~ 10 could have much longer lifetime 
in the phase of < 0.5 ms. However, for lower ellipticity, their life- 
times in the phase of < 0.5 ms are long enough for quark stars with 
~ IMq. Once a pulsar with spin period < 0.5 ms is ever found, 
low mass quark stars will be physically identified. 

3.4 Spin-down rate P for newborn quark stars and neutron 
stars 

We also use Eq. (20) to calculate the period derivative (P) for the 
nascent sub-millisecond quark stars and neutron stars. 

Figure 3 is a P — P diagram that shows the spin-down evo- 
lution for quark stars with different masses. It is found that, for 

2 |http : / /www| . atnf.csiro.au/research/pulsar 
/catalogue/ 
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Figure 2. The relation of lifetime (in the phase of < 1 ms) and gravita- 
tional ellipticity e e for quark stars with masses of 0.00 IMq (solid line), 
O.OIMq (dot-dash line), 0.1M Q (dashed line), magnetic field B = 10 8 G 
and the bag constant j3 = 60MeV fm -3 . The lifetime in the phase of 
sub-millisecond period is shorter if the quark star's mass is higher. 



different ellipticity there are different properties. For high elliptic- 
ity such as e e = 10 _J , the P can be changing about ten orders 
of magnitude for different periods (see the steep slopes of dot-dash 
lines and dashed lines). The rotational energy losses in this case 
are dominated by the gravitational wave (GW) radiation. For low 
ellipticity such as e e = 10 , in most cases, the rotational energy 
losses are dominated by magnetic dipole (EM) radiation and the P 
changes with periods relatively slow (solid lines). 

As a comparison, we also calculate the period derivative (P) 
of a neutron star (with an initial period 0.5 ms, mass of \AMq 
and radius of 10 6 km). The results are shown in Figure 4. One can 
see that the P is changing with periods as large as ten orders of 
magnitude. It is found that the neutron star spins down much more 
quickly than low mass quark stars, because of neutron star's high 
mass (~ Mq) for higher efficiency of GW radiation. 



4 SUB-MILLISECOND PULSARS FORMED THROUGH 
ACCRETION IN BINARY SYSTEMS 

There is also an important mechanism of "spin-up in binaries" for 
sub-millisecond pulsars' formation, which is widely discussed in 
the literatures. We regard this as "sub-Keplerian case" and make a 
comparison with our proposed AIC model "super-Keplerian case". 
In this section, we will find the minimal periods for both neutron 
stars and bare quark stars spun up by accretion in binary systems. 
We assume that the initial rotational periods of newborn pulsars 
could have an "equilibrium period" with two characteristic param- 
eters: magnetospheric radius and corotation radius. The magneto- 
spheric radius (r m ) is the radius where the ram pressure of particles 
is equal to the local magnetic pressure, i.e. 

^ ^ 2 mM 3/ \ 2 /7 .( BpR 6 2/7 

rm = * Ra = = H mmM ] 
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Table 3. The minimal equilibrium period for quark stars and lifetimes due to GW and EM radiation in the phase of sub-millisecond 
period for quark stars with different masses (1O" 3 M , O.IMq, 1.4M ) in the sub-Keplerian case. t\ , T2 , T3 are calculated by using 
e e = 10 -6 , while fi, T2, ts are calculated by using e e = 10 -9 . The bag constant f3 is in unit of Mev fm - 3 , the accretion ratio a 
is in unit of the Eddington accretion rate M^dd- 



p 


a 


Bo(10 8 G) 


-Pcqmm(ms) 


ti (yr) 


n(yr) 


T"2(yr) 


T"2(yr) 


t"3 (yr) 


T"3(yr) 


60 


0.71 


1.1 


0.613 


2.9 x 10 7 


2.8 x 10 10 


1.3 x 10 4 


4.1 x 10 9 


1.7 x 10 2 


1.5 x 10 8 


110 


0.85 


1.4 


0.453 


5.1 x 10 7 


3.6 x 10 10 


2.3 x 10 4 


5.6 x 10 9 


2.3 x 10 2 


2.6 x 10 s 



Table 4. The minimal equilibrium period and lifetimes for GW and EM radiation in the phase of sub-millisecond period of different 
EOSs of normal neutron stars in the sub-Keplerian case. The mass and radius data of neutron stars are obtained from Figure 2 of 
Lattimer et al. (2004). r and f are lifetimes within sub-millisecond period for neutron stars using e e = 10 — 6 and e e = 10 — 9 
respectively. 



EOS Poqmin(ms) Mass(M ) Radius(km) B (10 8 G) M (10 17 g s" 1 ) r(yr) f(yr) 



AP4 


0.55 


2.21 


10 


2 


6.36 


103 


1.02 X 10 8 


GS1 


0.52 


1.38 


8.27 


2 


5 


255 


2.37 X 10 s 


PAL6 


0.60 


1.48 


9.24 


2 


6.38 


177 


1.65 x 10 s 


MS0 


0.76 


2.76 


13.31 


1 


2.91 


35 


3.48 x 10 7 


GM3 


0.75 


1.56 


10.93 


2 


9.47 


94 


8.54 x 10 7 


MSI 


0.76 


1.81 


11.67 


1 


2.58 


90 


6.83 x 10 7 



3.24 x 10V 
1.857 x 10 6 < 



iU 7 m- 1/7 r1 2/7 m; 



2/7 
-2/7 



(22) 



where fj, m is the magnetic moment of per unit mass of the compact 
star; Bs is the surface magnetic strength in units of 10 8 G and Mn 
is the accretion rate in units of 10 17 g s _1 ; <j> ls me rat i° between 
the magnetospheric radius and the Alfven radius (Wang 1997; Bur- 
den & King 1998). Wang (1997) studied the torque exerted on an 
oblique rotator and pointed out that <j> decreased from 1.35 to 0.65 
as the inclination angle increased from 0° to 90°. Here we take 
<f> ~ 1, the influence of <f> is discussed in §6. 

When r m is very close to the compact star's radius, we could 
rewrite the accretion rate M in units of Eddington accretion rate 
(AfEdd), with a ratio, a, so that 



M — aMEdd = a 



= 1.0xl0 18 aA/ 1 / 3 /3- 1 / 3 gs- 1 .(23) 

With these equations obtained above, then we can get r m for quark 
stars, 



9.6a 



-2/7 R 4/7 J w-l/3 fl -10/21 



km. 



(24) 



S 1/3 2 / ^ 

The corotation radius is r c = 1.5 x 10 M t P ' cm. The spin 
periods of compact stars cannot exceed the Kepler limit via accre- 
tion. When the compact star was spun up to the Kepler limit by 
the accreted matter falling onto the compact star's surface, for neu- 
tron stars, as the equatorial radius expanded, one can use the simple 
empirical relation for the maximum spin frequency 



Onax = 7700Mi /2 R~ 3/2 s" 1 



(25) 



(Haensel & Zdunik 1989; Lattimer & Prakash 2004), which leads 
to 



P cq ^ 0.816M~ 1/2 .R 3/2 ms, 



(26) 



where M and R refer to the neutron star's mass and radius of non- 
rotating configurations. 

For quark stars, Gourgoulhon et al. (1999) used a highly 
precise numerical code for the 2-D calculations, and found that 
the f2 m ax could be expressed as f2 max = 9920\/^60 rad s , 



where /3 6 o = /V(60 MeV fm 3 ), which implied that P eq ^ 
0.633/3 14 1//2 ms. These are the so-called "sub-Keplerian condi- 
tion". 

The accretion torque, N, exerted on the compact star con- 
tains two contributions: one is positive material torque which is 
carried by the materials falling onto the star's surface; the other is 
magnetic torque which can be positive or negative, depending on 
the fastness parameter lu s = fi^/f^K = {r m /r c ) 3 ^ 2 . It is sug- 
gested that all the torques may cancel one another if the fastness 
is lj s = (r m /r c ) 3/2 w 0.884 (Dai & Li 2006). This implies a 
magnetospheric radius of r m = 0.92r c w r c . One can obtain an 
equilibrium period of P cq when setting r m = r c , 



0.512B 6 8 /7 fc 4 5/7 a- 3 ? 7 



(a) 



1 S17QBti 7 M^ /7 Rl 8/7 M^ /7 ms. (b) (2?) 



For quark stars, the equilibrium period is independent of mass 
and radius, and only dependent on bag constant, surface magnetic 
field, and accretion rate. Take Bo in the range [10 8 G, 10 12 G], 
we may use Eq. (27a) to calculate the minimal equilibrium pe- 
riod of different EOSs (equation of state) for quark stars. For 
j3 = 60 MeV fm" 3 , when a = 0.71, Bg = 1.1, one can get 
the minimal period 0.613 ms. For /3 = 110 MeV fm" 3 , when 
a = 0.85, Bg = 1.4, the minimal period is 0.453 ms. (See results 
in Table 3.) 

For neutron stars, data for mass and radius in different EOSs 
were taken from Lattimer & Prakash (2004, their Figure 2), So is 
in the range [10 8 G, 10 12 G]. The minimal equilibrium period is 
calculated using the Eq. (27b). (See results in Table 4.) 

In the sub-Keplerian case, the timescales in the phase of sub- 
millisecond for quark stars of different mass and neutron stars of 
different EOSs are listed in Table 3 and Table 4, respectively. For 
typical quark stars as well as neutron stars with high e e , their life- 
times in the phase of sub-millisecond period are about 10 2 years, 
which result in a too low detection possibility. However, for low e e , 
the lifetime of a sub-millisecond pulsar (even with a high mass) is 
long enough. 
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log 10 (period (sec)) 

Figure 3. Spin-down evolution of quark stars due to GW and EM radiations 
(period derivative versus spin period), with masses of O.IMq, 0.0 lMg, 
O.OOIMq. We choose ellipticity to be 10~ 5 (dot-dash lines), 10~ 7 (dashed 
lines), 10 — 9 (solid lines) in the calculation. It is evident that GW radiation 
dominates for quark stars with higher e e , while EM radiation dominates for 
lower e e . 



5 CONCLUSIONS AND DISCUSSIONS 

If a sub-millisecond pulsar is ever found, we have shown that it 
could be a quark star based upon plausible scenarios for its origin, 
the energy available for radiation and its lifetime. A new possible 
way to form sub-millisecond pulsars (quark stars) via AIC of white 
dwarfs has been discussed in this paper. In the super-Keplerian 
case, we derived the initial period P q via angular momentum con- 
servation with consideration of the special and general relativistic 
effects, and calculated the lifetime and gap parameters of a new- 
born quark star. Quark stars with different masses could have the 
minimal rotational period around 0. 1 ms. In most cases, quark stars 
would be bare (Xu 2002), therefore, a vacuum gap would be formed 
in the polar cap region. Based on our rough estimations without 
considering the effect of frame dragging (Harding & Muslimov 
1998), we found that the basic parameters (including rotational en- 
ergy loss) in the gap are suitable for pair (electrons and positrons) 




Figure 4. Period derivative versus spin period diagram for a neutron star 
with an initial period of 0.5ms, mass of 1.4M and radius of 10 6 km. The 
neutron star spins down quickly due to high mass (moment of inertia) for 
GW radiation. 



production and sparking. They can be detected as sub-millisecond 
radio pulsars. 

We also used an approximate formula to calculate nascent 
quark star's moment of inertia, but there are no accurate solutions to 
fast rotating compact stars' configuration until nowadays. It should 
be investigated precisely in the future. In the calculation of WD's 
mass and radius, we just considered the non-rotating configuration. 
But it does not change the conclusions of this paper. If the central 
density p c of the WD is lower than 10 11 g cm~ 3 before collapsing, 
the resulting WD has a larger radius and moment of inertia, con- 
sequently, the newborn quark star could have a smaller spin period 
«lms). 

Both the special and general relativistic effects are weak for a 
low mass (e.g. Jupiter-like) quark star with a small radius. The ro- 
tational energy is lost via GW and EM radiation. The GW radiation 
dominates the rotational energy loss in the phase of sub-millisecond 
period, if magnetic field of stars is not so large. Such quark 
stars therefore have long lifetimes (several million years if mass 
~ lO~ 3 Af0) to maintain their spin periods of sub-millisecond. We 
have considered the bar-mode of GW radiation in this paper, while 
other GW mo de (e.g. r-m ode) may be important but not yet con- 
sidered here JXuH2006bh . The subsequent relaxation timescale of 
a newborn quark star to a rigidly rotating configuration could be 
negligible since a quark star may be solidified soon after birth. 

An important constraint for sub-millisecond pulsar's detection 
is its lifetime in the phase of < 1 ms due to GW and EM radiation. 
A possible method is proposed to constrain the lower limit of the 
pulsars' equatorial ellipticity, i. e., e , m i n ~ 10~ 9 , by evaluating 
millisecond pulsars' period derivative via Eq. (20). For larger el- 
lipticity, e. g., £ e = 10 -6 , it is clear that GW radiation dominates 
the energy loss in the phase of short period for either recycled or 
AIC's compact stars. The corresponding lifetime is shorter for a 
compact star with higher mass (~ Mq), but longer for a star with 
lower mass (~ O.OOIMq). However, if the ellipticity is lower, e. 
g., e e = 10 , EM radiation dominates the rotational energy loss. 
The corresponding lifetime of a quark star (even with a high mass 
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~ Mq) is long enough, and there are no lifetime constraints for 
sub-millisecond pulsars' detection. Solid evidence of quark stars 
will be obtained if a pulsar with a period of less than ~ 0.5 ms is 
discovered in the future. 

In the sub-Keplerian case, neutron and "bare" quark stars can 
be spun up to sub-millisecond periods (even ~ 0.5 ms) through 
accretion in binary systems. When neutron stars are spun up to 
the Kepler limit, the minimal equilibrium periods depend only on 
the mass and radius of the nonrotating configurations. Quark stars' 
minimal equilibrium periods depend on the bag constant. 
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APPENDIX A: DIFFERENTIALLY ROTATING WD 
MODEL 

WD could be rotating differentially. As stated by (Mueller & 
Eriguchi 1985), the WD's angular velocity Q is a function of the 
distance from the rotation axis ui. The angular momentum distribu- 
tion (so-called rotation law) is 

W^J&FTW' (A1) 

where Q. c is the central angular velocity, R e is the equatorial radius, 
and a is a free parameter. When differential rotation is taking into 
account, we can numerically evaluate the angular momentum of the 
WD's inner collapsed core, i.e., 



Jcoro = /J Ji = 2~] / cr^nrf sin 3 9Q(ri sin 
i i J o 



(A2) 



i + 



(r t -0.5, 2 ' — g a i?wDln 



,, a 2 R 2 Tr 1 
2 n 2 V WD ' i 



1 - 



where Ji is the angular momentum of each spherical shell with 
integral radius n. 



